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We consider the possibilities for obtaining information about the equation of state for quark 
matter by using future direct observational data on gravitational waves. We study the nonradial 
oscillations of both fluid and spacetime modes of pure quark stars. If we observe the / and the 
lowest wii modes from quark stars, by using the simultaneously obtained radiation radius we can 
constrain the bag constant B with reasonable accuracy, independently of the s quark mass. 
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^ . I. INTRODUCTION 

D 

' Towards observing gravitational waves, several gravitational wave interferometers on Earth are steadily developing 
today, such as the Laser Interferometric Gravitational Wave Observatory (LIGO) Q, TAMA300 0, GEO600 and 
VIRGO Q . Thus it will be possible for us to detect gravitational waves directly in the near future. It is believed that 
mergers of binary neutron-star-neutron-star (NS-NS), neutron-star-black-hole (NS-BH), and BH-BH, or supernovae, 
04 and so on, can become strong sources of gravitational waves. After these violent events occur, compact objects may be 
^ ' left and may be turbulent. Then gravitational waves are emitted from them. At this time, these gravitational waves 
, convey information on the source object. If these gravitational waves are directly detected on Earth, it is possible to 

■ obtain some information about the sources. This research field is called "gravitational wave astronomy." In this field, 
there is an attempt to obtain information about properties of the equation of state (EOS) of high density matter. 

I This is one of the most important purposes of gravitational wave astronomy. 
\ Gravitational waves are emitted by the nonspherical oscillation of compact objects. The oscillations are damped out 

■ as gravitational waves carry away the oscillational energy. Such oscillations are called quasinormal modes (QNMs). 
I The QNMs have complex frequencies whose real and imaginary parts correspond to the oscillational frequency and 

damping rate, respectively. The QNMs fall into two types from their nature. One involves fluid modes which are 
^ ' connected with the stellar matter. The other involves spacetime modes which are the oscillations of spacetime metric. 
f^j^, Moreover, the fluid modes are classified into various types. The well-known modes are the /, p, and g modes 

The fluid modes have a characteristic that the damping rate Im(w) is much smaller than the oscillational frequency 
Re(a;). The / mode is the fundamental mode. There exists only one / mode for each index I of spherical harmonics 
Yim- The p mode is the pressure or acoustic mode, whose restoring force is caused by the pressure gradient inside 
^ I the star. The g mode is the gravity mode, which arises from buoyancy in a gravity field. The w and wu modes are 
spacetime modes 0,0- Unlike the fiuid modes, the damping rate of the w and wu modes is comparable to or larger 
than the oscillational frequency. 

For the QNMs of neutron stars, so far many authors have argued the possibilities for determining the EOS in 
the high density region and/or for restricting the properties of neutron stars, such as the radius R or mass M, by 
employing the observed gravitational wave of several nonradial modes |E B 0. ^> Oi Qi ■ As a candidate for 
a star which is smaller than neutron stars, the possibility of a quark star or a compact star, which is supported by 
degenerate pressure of quar k matter, has been pointed out. Such a quark star has been investigated by many authors 
(see, e.g., Refs. [H El US lIlllilMilllllMlillSIMIlllMllSllSllllll references therein). In their 
view, it is commonly assumed that such quark stars contain quark matter in the core region and are surrounded by 
hadronic matter, although they are in the branch of neutron stars '28'|. Witten suggested another type of quark 
star. If the true ground state of hadrons is bulk quark matter, which consists of approximately equal numbers of u, 
d, and s quarks ("strange matter"), there exist self-bound quark stars. They are called "strange stars." In this case, 
their mass and radius are smaller than those of typical neutron stars, which are ~ 10 km and ~ 1.4Mq, respectively. 
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Because we do not have reliable information about the equilibrium properties of hadronic and quark matters at high 
densities, it is not clear what kinds of quark stars are realized. 

Recently, Drake et al. reported that the deep Chandra LETG+HRC-S observations of the soft X-ray source RX 
J1856. 5-3754 reveal an X-ray spectrum quite close to that of a blackbody of temperature T = 61.2 ± 1.0 eV 5^]. The 
data contain evidence for the lack of spectral features or pulsation |3^. Drake et al. also reported that the interstellar 
medium neutral hydrogen column density is Nh = (0.8-1.1) x 10^° cm~^. With the results of recent HST parallax 
analyses, that yields an estimate of 111-170 pc for distance D to RX J1856. 5-3754. Combining this range of D with 
the blackbody fit leads to a radiation radius of i?oo = 3.8-8.2 km. That is smaller than typical neutron star radii 3J^. 
Thus they suggested that the X-ray source may be a quark star. 

In the meanwhile, in Ref. [s^, Walter and Lattimer claimed that the blackbody model adopted in Ref. ||33] could 
not explain the observed UV-optical spectrum. They undertook to fit the two-temperature blackbody and heavy- 
element atmosphere models which were discussed in Ref. in X-ray and UV-optical wavelengths. From their 
analyses they found that the radiation radius was 12-26 km and was consistent with that of a neutron star. However, 
this model cannot explain the lack of spectral features. In addition, recently Braje and Romani also suggested a 
two-temperature blackbody model, which can reproduce both the X-ray and optical- UV spectral 38]. Although their 
model is inconsistent with the fact that pulsation is not detected, this might be explained by the object being a young 
normal pulsar and its nonthermal radio beam missing the Earth's line of sight. However, they can not also answer 
why there are no features in the observed X-ray spectrum. In Ref. ^S^l, Burwitz et al. discussed the possibility of a 
condensate surface which is made of unknown material to explain both the UV-optical and X-ray spectra in a neutron 
star. More recently, some groups discussed the possibility that the effects of a strong magnetic field {B ^ 10^^ G) ^3 
or rapid rotation (P < 10 ms) ^^^Y smear out any spectral features. In these situations, however, it seems that 
there are no reliable models that account for all the observational facts. It is still controversial whether RXJ1856.5- 
3754 is a normal neutron star or some other compact star like a quark star. Here we adopt the simplest picture, which 
is used by Drake et al. |33l |: the uniform temperature blackbody model. 

As for gravitational waves emitted from quark stars. Yip, Chu, and Leung studied nonradial stellar oscillation for 
stars whose radius is around 10 km |42j |. Kojima and Sakata demonstrated the possibility of distinguishing quark 
stars from neutron stars by using both the oscillational frequency and the damping rate of the / mode ,43] . Sotani 
and Harada showed that the / and lowest wu modes depend strongly on the EOS of quark matter and the properties 
of quark stars, where the lowest wu mode is the one which has the largest frequency among all wu modes j44j . Then 
they pointed out the possibility of determining the EOS and/or the stellar properties. Furthermore, they also studied 
w modes in detail. However, they clearly showed that w modes do not depend much on the EOS of quark matter and 
are not important for constraining the model parameters from observations. In their work, however, they assumed 
that the star is a pure quark star and that the EOS is described by a simple bag model which has only one parameter: 
i.e., the bag constant B. In general, there are a variety of parameters even within bag models: e.g., the bag constant 
B, the strange quark mass rus, the fine structure constant in QCD ac, and so on. In particular, if the effects of 
nonvanishing strange quark mass rUs are taken into account, the structure of quark stars can be affected considerably 
(for recent analyses, see Ref. IZsj and references therein). In this situation, we compute the QNMs in the bag model 
used in Ref. citeKohri:2002hf and investigate the possibility of restricting the model parameters by the observations 
of QNMs. In this study, we deal with only / and the lowest wu modes in response to the results in Ref. [3 ]. 

Effective theories of quantum chromodynamics (QCD), such as bag models, perturbation theories, or finite- 
temperature lattice data, are difficult to test by experiments, especially in a low temperature and high density 
regime. To further study them and fit their model parameters, we should compare the theoretical predictions in such 
models with experimental data — e.g., data in relativistic heavy ion collision experiments and so on. Thus, in this 
situation information about compact objects obtained by astrophysical observations is indispensably valuable, being 
independent of the above-mentioned ground-based experiments in nuclear physics or particle physics. Among the 
astrophysical observations, the observation of gravitational waves emitted from oscillating compact objects is quite 
unique because gravitational waves directly convey information on the internal structure of compact objects, where 
the density would reach the nuclear density. 

The plan of this paper is as follows. In Sec. ^we introduce the basic equation including the EOS to construct 
quark stars as the source of gravitational waves and to give the properties of quark star structure. We present a 
method of determining the QNMs for the case of spherically symmetric stars in Sec. IIIII In Sec. IIVI we show the 
numerical results for the QNMs for the quark star constructed in Sec. ^1 In this section we present the dependence 
of the QNMs on the parameters of the EOS and stellar properties, and discuss the possibility of determining the EOS 
of quark matter. We conclude this paper in Sec. O We adopt units oi c = h = G = 1, where c, h, and G denote 
the speed of light, reduced Planck's constant, and gravitational constant, respectively, and the metric signature of 
(— , +, +) throughout this paper. 
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II. QUARK STAR MODELS 

We assume that the quark star is static and spherically symmetric. In this case the metric is described by 

ds^ = -e^'^dt^ + e^^dr^ + {d0^ + sin^ Odct)^) , (2.1) 
where $ and A are metric functions of r and A is related to the mass function ■m{r) as 

m(r) ^^r{l- e-^A) . (2.2) 

The mass function m{r) is the gravitational mass inside a surface of radius r and satisfies 

dm ,9 „, 
— (2.3) 

where p is the energy density. The equilibrium stellar model is constructed by solving the equations 

dP [p + P){m + inr^P) 



dr r (r — 2m) 

d$ (m + Airr^P) 
dr r (r — 2m) 



(2.4) 
(2.5) 



Equation H2.4|) is the Tolman-Oppenheimer-Volkoff (TOV) equation In addition to the above equations we need 
an EOS to compute the star configuration. The stellar surface r = i? is the position where the pressure vanishes, and 
the stellar mass M is defined as M — m{R). Because the metric in the exterior of the star is the Schwarzschild one, 
the metric (12. 1|) must be connected smoothly with the Schwarzschild one at the stellar surface. 

Next we describe the EOS in bag models. To construct stars composed of zero-temperature uds quark matter, we 
start with the thermodynamic potentials fi^ for a homogeneous gas of q quarks {q = u, d, s) of rest mass ruq up to 
first order in the QCD coupling constant ac- The expressions for fL are given as a sum of the kinetic term and the 
one-gluon-exchange term at the renormalization scale A = niq j24ll47ll48l |. Here we assume that m„ = to^ = 0. Then 
ilq is expressed by 



47r2 



TT 



47r2 



TT 



m'^ ( 3 3 2a 



TT 



3F{xs) {F{xs) + 2 In a;,) - 2??4 + 6 ( In A ) f{xs) 



(2.8) 



with F{xs) = XsTjs — In (xs +'^s)^ where Xs — Ps/mg, rjs = A/a;^ — 1, pg is the chemical potential of q quarks, 
Q!c = Off, (A), and rrig — mg{A). Here we choose A = pg. The electron thermodynamic potential in the massless 
noninteracting form is given by 

n. = (2^9) 

where pe is the electron chemical potential. 

Within the bag model, we express the total energy density p as 

p= {n, + p,n,) + B, (2.10) 

i—u,d,s,e 

where B is the bag constant |52l | which is the excess energy density effectively representing the nonperturbative color 
confining interactions, n.i is the number density of i particles given by rii = —dQi/dpi. Then, the pressure p is 
represented by 

p = - J2 ^—B. (2.11) 

i—u,d,s,e 
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The parameters nis, ac, and B were obtained by fits to light-hadron spectra; e.g., see Refs. |49l. l5Ct l5l| (also see 
Table nj. Note that these parameters do not correspond directly to the values appropriate to bulk quark matter [s^ . 
When we allow for possible uncertainties, we set 145 < B'^^'^ < 200 MeV, < < 300 MeV, and < ac < 0.9. The 
ranges of nis and ac are not always consistent with the values in Table However, when we renormalize them at an 
energy scale of interest to us, they are consistent with values of the Particle Data Group js^. 

To obtain the equilibrium composition of the ground-state matter at a given baryon density or chemical potential, 
the conditions for equilibrium with respect to the weak interaction and for overall charge neutrality are required. 
These conditions are expressed by 

/i« = ^J■d - Me, (2-12) 
= Ms, (2-13) 

and 

2 11 

-^riu - -Ud - -Us = He- (2.14) 

From Eqs. (|2. 6(1 - (|2. 14(1 . we can obtain the energy density and pressure as a function of ub (or M-b). In Ref. 0| the 
authors assumed that mg ~ and ac ~ 0. Then, the EOS is analytically given by 

p = \{p-m- (2.15) 

In this case, the density at the stellar surface is given by p(K) — AB. 

In Fig. ^ we plot several sequences of quark star models with different values of the bag constant B and the s quark 
mass nis- Here we have fixed the fine structure constant in QCD as ac = 0.6 because the dependence of properties 
of quark stars on ac is weak, as indicated in Ref. |4^. Therefore we deal with stellar models only for the case of 
ac — 0.6 in this paper. In Fig. QJa) we plot the mass M/Mq of the star as a function of central density pc- Here we 
set B^^^ = 145, 200 MeV and nis = 0, 150, 300 McV, respectively. In Fig. Ub) we also plot the radius of the star 
as a function of central density by thick lines. Thin lines denote the radiation radius Roo- In Fig. ^ the solid line, 
dashed hue, and dotted hne correspond to the cases of = 0, 150, and 300 MeV, respectively. In Tables HTl and IIIII 
we list several properties of quark stars whose radiation radii are within the range of Roc = 3.8-8.2 km, which was 
reported by Ref. .33.] . 



III. DETERMINATION OF THE QNM 

Since we are interested in the dependence of QNMs on the EOS of quark matter, we consider only polar perturbation. 
If we adopt the Regge- Wheeler gauge, the metric perturbation is given by 



Or 



,2A 








V 



r'+^i^sin^i 



^ Tin ^ , 



(3.1) 



where gj^"^ is the background metric 1(2.1(1 of a spherically symmetric star and H, Hi, and K are perturbed metric 
functions with respect to r. We apply a formalism developed by Lindblom and Dctweiler "54] for relativistic nonradial 
stellar oscillations. The components of the Lagrangian displacement of fluid perturbations are expanded as 



= - 



sin^ I 



(3.2) 
(3.3) 
(3.4) 



where W and V are functions of r. 
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Assuming that the matter is a perfect fluid and the perturbation is adiabatic, we have the following perturbation 
equations derived from Einstein equations: 



dHi 
dr 



dK 
dr 

dW 

dr 

dX_ 
dr 



2m OA 
; + 1 + e^^ 

r 



47rr^(P - p)e 



2A 



,2A 



r 

2r 

l + l 



H + k + 16tt{P + p)V 



Hi + -H- 
r 



l + l d$ 



W 



1 



I 



r^ 2 
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[^dr r y 





^(^ + 1) 

2r 



Hi 



1 

2^ 



2 



K 



l{l + l)d^ 



dr 



V 



d 
dr 



dr J 



W 



3 m 
r 



1(1 + 1) 



+r'^e 



2 

-2A 



2 -2* 

w e 



1 I 2 2 -2* 

1 + uj r e 



l{l + l)d<P 
2r dr 
l{l + l) 



- r - 3m - A'nrT) — 
a?- 



o2* 



V 



uj^{P + p) 



1 dP 



1 



(P + p)H 



r dr 2 ' 

where 7 is the adiabatic index of the unperturbed stellar model, which is given by 

p + PdP 



1 ■ 



P dp ' 



(3.5) 
(3.6) 

(3.7) 



(3.8) 



(3.9) 
(3.10) 

(3.11) 



Furthermore, if the proper boundary conditions are imposed on the above perturbation equations, the problem to 
solve becomes an eigenvalue one. These boundary conditions are given as follows: (i) the eigenfunctions are regular at 
the stellar center, (ii) the Lagrangian perturbation of pressure vanishes at the stellar surface, and (iii) the gravitational 
wave is only an outgoing one at infinity. The QNMs of quark stars are determined by solving this eigenvalue problem, 
and UJ is the eigenvalue of the above perturbation equations. For the treatment of the boundary condition at infinity, 
we adopt the method of continued- fraction expansion proposed by Leaver [ssj . The detailed method of determining 
QNMs is given in [illHfill. 



IV. NUMERICAL RESULTS 



A. / mode 

As for gravitational wave radiation emitted from a compact star, there exist QNMs which correspond to each I. 
In this paper, for simplicity we discuss such QNMs only for I — 2, which would be a dominant mode in a weak field 
regime. 

Compared with the other QNMs, we might relatively easily detect the / mode because its frequency and damping 
rate are approximately within the sensitivity band of gravitational wave interferometers: e.g., LIGO and so on. 
Therefore, first we discuss the possibility of constraining parameters of bag models by using the / mode data. When 
we discuss the observational data of frequency Re(a;) and damping rate Im(u;) of the / mode, we should bear in mind 
that the observational data of frequency will be far more accurate than that of damping rate. It was discussed in 
Ref. ^3 that the relative error of the frequency will be about three orders of magnitude smaller than that of the 
damping rate. To get an insight into the dependence of the / mode on the bag constant B and the s quark mass m^, 
we compute the / mode for each stellar model in Tables ITU and IIIll 
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In Fig. Ola) we plot the complex frequencies of / mode for 1 = 2. In the figure, the squares, circles, and triangles 
represent nis = 0, 150, and 300 MeV, respectively, and the sohd (open) marks denote the case of B^/^ = 145 MeV 
(200 MeV). In Ref. |4J| the authors did not consider such dependence because they simply assumed that iris = 
as their first step in the paper. From Fig. [2Ia), however, we find that both the frequencies and damping rates change 
at about 10% as nis changes from zero to 300 MeV for each B, and we see that it is important for us to include such 
effects of nonzero values of in bag models to predict the / mode. 

In Fig. I^lwe plot the radiation radius as a function of the (a) frequency Re(w) and (b) damping rate Im(a;) of the 
/ mode for 1 = 2. In each panel, the solid and open marks denote i?^/* = 145 MeV and 200 MeV, respectively. From 
Fig. Ol^b) we see that the damping rate is very sensitive to the radiation radius, unlike the frequency in Fig. |2Ia). 
If we observe the damping rate of the / mode within 20%, the radiation radius can be determined within 10% as a 
function of B, independently of the value of m^. On the other hand, from Fig. ISl^a) we see that if we observe the 
frequency within about 20%, for example, we can determine the value of B with reasonable accuracy by using the 
simultaneously obtained radiation radius from Fig. OJb). To demonstrate that, we investigate the dependence of the 
/ mode on both the bag constant and the s quark mass when we fix the radiation radius as Roo = 6.0 km. In Fig. 01 
we plot the contours of the (a) frequency Re(cj) and (b) damping rate Im((jj) of the / mode for ^ = 2 in the {B^^^, 
nis) plane. The numbers on the curves denote the value of frequency Re(ci;) or damping rate Im(a;) in units of kHz. 
We should note that each contour line is nearly parallel to the vertical axis in both Figs. Ela) and^^b). This implies 
that we can obtain a stringent constraint on the bag constant. Indeed, Fig. Efa) shows that if the frequency of the / 
mode is detected within 20%, then we can determine the value of the bag constant within about 15% independently 
of the value of the s quark mass. 

B. The lowest wu mode 

In general, there exist infinite number of w modes and several wu modes for each I. As was clearly shown by 
numerical analysis in Ref. j44| . however, the w modes are insensitive to the parameters in bag models in which we 
are interested here. Therefore, we do not adopt w modes in this study. In addition, only one wu mode was found in 
each quark star model adopted in this study. That is because the "compactness"- i.e., mass to radius ratio — of the 
adopted quark stars is much smaller than that of standard neutron stars (see Tables Hll and IIIl|) . Here, for simplicity 
we deal with only the lowest wu mode for I = 2 • 

In Fig. Olb) we plot the complex frequencies of the the lowest wu mode. From Fig. Hfb), we see that the nis 
dependence on the complex frequencies is important, although it was not considered in Ref. |4J|. In Fig. |Slwe plot 
the radiation radius as a function of the (a) frequency Re(cij) and (b) damping rate Im(a;). The notes on the marks 
are the same as those in Fig. 13 This figure tells us that it has the same tendency which appears in the case of the 
/ mode; namely, the damping rate is very sensitive to the radiation radius, unlike the frequency. Thus, if we directly 
observe the damping rate, we can constrain the radiation radius of the quark star. 

Furthermore, along with the case of the / mode, we plot the contour of the frequency and the damping rate of the 
lowest Wu mode in the {B^^^, rUs) plane in Fig. |H1 By using the lowest wu mode, we are able to develop a similar 
argument as that discussed in the case of the / mode and get an independent constraint on the bag constant and s 
quark mass. 

However, we should keep in mind that even if a large amount of energy is released through these modes, accurate 
observation of the w and wu modes may be difficult. That is because both the frequency and damping rate of these 
modes are larger than the sensitivity ranges of gravitational wave interferometers §j . Thus we would mainly use the 
observational data of the / mode to obtain information about the bag model parameters. Then we could subsidiarily 
use the data on the lowest wu mode. It should be noted that, since Fig. |S|has similar features found in the / mode, 
even if we combine it with information from the / mode, we may not have strict constraints on both TOs and -B, 
independently. There exists an essential degeneracy in / and the lowest wu mode QNMs because the configuration 
of the contour of Fig. 0]is very similar to that of Fig. 

V. CONCLUSION 

We have discussed how we can obtain information about the EOS of quark matter by using future observations of 
gravitational waves emitted from quark stars. In particular we have studied the EOS in bag models and assumed that 
the star is a pure quark star. We have computed the QNMs — i.e., the / and the lowest wu modes — in several quark 
star models. We have demonstrated that by comparing the results of theoretical computations with the observational 
data of the / mode and the lowest wu modes we can obtain constraints on the bag constant B and s quark mass m^. 
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If we have the damping rate — i.e., Im(a;) — of the / mode within 20%, the radiation radius of the quark star can 
be determined within about 10% including the uncertainty of the s quark mass. Furthermore if we also obtain the 
frequency — i.e., Re(a;) — of the / mode within 20%, the value of the bag constant can be determined within about 
15%, independently of the uncertainty of the s quark mass, by using the simultaneously obtained radiation radius. 
Concerning the lowest wu mode, we can also develop a similar argument as in the case of the / mode and get 
independent constraints on the model parameters. However, note that it is relatively difficult to detect the lowest 
Wu mode by the future planned gravitational wave interferometers whose frequency ranges are not very sensitive to 
the lowest wu mode. Therefore, such data will be subsidiarily used in statistical analyses. It should be also noted 
that there is a degeneracy in the dependence of / and the lowest wu mode QNM complex frequencies on the bag 
model parameters and B. As for high frequency gravitational waves, a dual-type detector has been proposed, 
which would reach very good spectral strain sensitivities (~ 2 x 10~^^ Hz) in a considerable broadband (between 1 
and 3 kHz) and then open a new interesting window to the QNMs of compact objects. 

As long as we have data with tolerable accuracy in future, we will be able to perform statistical analyses to fit 
them and get constraints on both the bag constant and the s quark mass. If we have further information about 
the bag constant or the s quark mass by utilizing the EOS models based on future developments in an effective 
theory of QCD, such as perturbation theories or finite-temperature lattice data, we can constrain them more strictly. 
Including independent observations of the radiation radius — e.g., sorts of X-ray observations — would also support us 
in inferring the bag model parameters. QNMs from compact stars will be observed in the near future and deepen our 
understanding of hadron physics and QCD. 
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TABLE L Bag model parameters fitted to hadron mass spectra. 



B^/* (MeV) 


(MeV) 


etc 


Reference 


145 


279 


2.2 


m 


200-220 


288 


0.8-0.9 


[50] 


149 


283 


2.0 


f51] 
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TABLE II: Properties of quark stars in the case of B^^"^ — 145 MeV. Here, -Roc, rns, Pc, M/Mq, R, and M/R are the radiation 
radius, the s quark mass, the central density, the mass in units of Mq, the radius, and the compactness of the star, respectively. 



Roo (km) nis (MeV) (g/cm^) M/Mq R (km) M/R 



3.800 


0.0 


4.221 X 10" 


4.539 X 10" 


2 


3.731 


1.796 X 10" 


2 


3.800 


150.0 


4.714 X 10" 


5.020 X 10" 


2 


3.724 


1.991 X 10" 


2 


3.800 


300.0 


5.027 X 10" 


5.305 X 10" 


2 


3.719 


2.106 X 10" 


2 


6.000 


0.0 


4.412 X 10" 


1.678 X 10" 


1 


5.735 


4.320 X 10" 


2 


6.000 


150.0 


4.979 X 10" 


1.846 X 10" 


1 


5.706 


4.778 X 10" 


2 


6.000 


300.0 


5.392 X 10" 


1.951 X 10" 


1 


5.688 


5.066 X 10" 


2 


8.200 


0.0 


4.731 X 10" 


3.925 X 10" 


1 


7.544 


7.684 X 10" 


2 


8.200 


150.0 


5.436 X 10" 


4.293 X 10" 


1 


7.472 


8.483 X 10" 


2 


8.200 


300.0 


6.049 X 10" 


4.540 X 10" 


1 


7.423 


9.031 X 10" 


2 



TABLE III: Same as Table El except for B^/"* = 200 MeV. 



fioo (km) nis (MeV) pc (g/cm'^) M/Mq R (km) M/R 



3.800 


0.0 


1.660 X 10^^ 


1.477 X 10" 


1 


3.560 


6.125 X 10" 


2 


3.800 


150.0 


1.796 X 10^^ 


1.563 X 10" 


1 


3.544 


6.513 X 10" 


2 


3.800 


300.0 


2.019 X 10^^ 


1.692 X 10" 


1 


3.520 


7.098 X 10" 


2 


6.000 


0.0 


2.038 X 10^^ 


4.742 X 10" 


1 


5.113 


1.369 X 10" 


1 


6.000 


150.0 


2.285 X 10^^ 


4.977 X 10" 


1 


5.052 


1.455 X 10" 


1 


6.000 


300.0 


2.781 X 10^^^ 


5.344 X 10" 


1 


4.953 


1.593 X 10" 


1 




(a) (b) 

FIG. 1: (a) Mass of the star as a function of central density and (b) radius of the star as a function of central density (thick 
lines) for B^^"^ = 145 MeV and 200 MeV. Thin lines denote the radiation radius -Roo. Here we adopt three cases of the s quark 
mass: i.e., rus = (solid line), 150 MeV (dashed line), and 300 MeV (dotted line). We adopt Qc = 0.6 as a value of the fine 
structure constant in QCD. 
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FIG. 2: Complex frequencies of (a) / mode and (b) the lowest wu mode for Z = 2. The squares, circles and triangles represent 
rus — 0, 150 MeV and 300 MeV, respectively. In each panel, the solid (open) marks denote the case of B^^'^ = 145 MeV (200 
MeV). 
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FIG. 3: Radiation radius as a function of (a) frequency Re(a;) and (b) damping rate Im(ii;) of / mode for 1 = 2. In each figure, 
the solid and open marks denote B^^* = 145 MeV and 200 MeV, respectively. 
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FIG. 4: (a) Contours of frequency Re(aj) and (b) damping rate Im(u;) of / mode for Z = 2 in {B^'''^, rus) plane. Here we adopt 
the radiation radius as Roa ~ 6.0 km. The numbers with the curve denote the value of frequency Re(ci;) or damping rate Im(ci;) 
in units of kHz. 
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FIG. 5: Radiation radius as a function of (a) frequency Re(ti;) and (b) damping rate Im(a') of the lowest wu mode for I = 2. 
In each figure, the solid marks and open ones denote B^^"^ = 145 MeV and 200 MeV, respectively. 
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FIG. 6: (a) Contours of frequency Ke{iij) and (b) damping rate Im(a;) of the lowest wu mode for Z = 2 in (-B^^*, rus) plane. 
Here we adopt the radiation radius as Roo = 6.0 km. The numbers with the curve denote the value of frequency Re(ti;) or 
damping rate Im(aj) in units of kHz. 



